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A direct classical analog of the quantum dynamics of intrinsic decoherence in Hamiltonian sys- 
tems, characterized by the time dependence of the linear entropy of the reduced density operator, 
is introduced. The similarities and differences between the classical and quantum decoherence dy- 
C^) . namics of an initial quantum state are exposed using both analytical and computational results. 

In particular, the classicality of early-time intrinsic decoherence dynamics is explored analytically 
using a second-order perturbative treatment, and an interesting connection between decoherence 
rates and the stability nature of classical trajectories is revealed in a simple approximate classical 
theory of intrinsic decoherence dynamics. The results offer new insights into decoherence, dynamics 
' of quantum entanglement, and quantum chaos. 
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INTRODUCTION 



Quantum dynamics induces unitary transformations in Hilbert space, but most often it is only the dynamics 



projected onto a Hilbert subspace that is of interest. In general this reduced dynamics is nonunitary and therefore 
displays decoherence pj. For example, if a system of interest is coupled to a bath, then averaging over the bath degrees 
\Q ' of freedom introduces decoherence in the system dynamics. Likewise, in an isolated system, the reduced dynamics of 
. a subsystem of this isolated system can display decoherence. We have termed decoherence in the latter case "intrinsic 
C*~) ' decoherence" since it does not involve an external bath 0- 

Understanding decoherence is of crucial importance to a variety of modern fields such as quantum information pro- 
cessing |3| and quantum control of atomic and molecular processes |4lae1. Our interest here is in the quantum-classical 
correspondence (QCC) between classical and quantum descriptions of the dynamics of decoherence. Specifically, we 
+Jj ■ consider an initial quantum state subjected to either quantum or classical dynamics and compare the time evolution 
of the decoherence in both cases. We note that the formal theory of correspondence between quantum dynamics 
and classical Liouville dynamics suggests that classical Liouville dynamics projected onto a subspace should also 
display decoherence. That is, as in the quantum case, the classical Liouville dynamics considered in the entire phase 
space is unitary and the classical Liouville dynamics projected onto a subspace is nonunitary. We therefore expect 
that the reduced classical Liouville dynamics propagated classically will show decoherence dynamics that is, at least 
qualitatively, parallel to that seen in the reduced quantum dynamics insofar as the loss of phase information, entropy 
production, etc. In the case of bath-induced decoherence we recently showed analytically that (a) one can indeed 
introduce a direct classical analog of quantum decoherence, and (b) examining the dynamics of decoherence classically 
gives new insights into both the dynamics of decoherence described quantum mechanically and into the conditions 
for QCC in decoherence dynamics 

Here we extend these considerations to intrinsic decoherence, both analytically and computationally. Specifically, in 
this paper we study QCC in the dynamics of intrinsic decoherence in smooth Hamiltonian systems, with an emphasis 
on the usefulness of classical dynamics in describing intrinsic decoherence. In particular, the classicality of early-time 
intrinsic decoherence dynamics is studied using a second-order perturbative treatment, and the interesting connection 
between decoherence rates at later times and the stability properties of classical trajectories is revealed by considering 
a simple approximate classical theory of intrinsic decoherence dynamics. The analytic and computational results shed 
new light on decoherence, dynamics of quantum entanglement, and quantum chaos. This study is also of interest 
to semiclassical decoherence studies Q, e.g., semiclassical descriptions of intrinsic decoherence dynamics in large 
molecular systems Q. 

This paper is organized as follows. In Sec. ^ we introduce a second-order perturbation theory in an effort to 
understand QCC in early-time intrinsic decoherence dynamics. For simplicity we focus upon two degree-of-freedom 
systems, but the extension to larger systems is straightforward. Computational results of two sample cases in coupled- 
oscillator model systems, which strongly support the physical picture afforded by the perturbative treatment, are 
presented in Sec. IIIII Then, a classical theory of intrinsic decoherence dynamics for initially localized states is derived 
in Sec. IIVI In the same section, detailed computational studies using this simple theory are carried out for the quartic 
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oscillator model and one of its variants. Discussions and a summary comprise Sec. [V] 

II. EARLY-TIME INTRINSIC DECOHERENCE DYNAMICS 

Consider a conservative system composed of two subsystems, with the total Hamiltonian given by 

H(Q, P, q,p) = ^- + V - + Vi(Q) + V 2 (q) + V 12 (Q, q), (1) 

where (Q,P) and (q,p) are dimensionless phase space conjugate variables, Vi is the potential of the i-th subsystem, 
and V\2 (Q,q) describes arbitrary coupling between the two subsystems. As the system evolves the total system 
wavefunction \ip(t)} becomes inseparable due to quantum entanglement, even if it is initially separable in Q and q. 
As a result, measuring a subsystem would collapse the system wavefunction and therefore affect the properties of the 
other subsystem. Similarly, ignoring a subsystem decoheres the other one. The degree of intrinsic decoherence, which 
is induced by, and is a manifestation of, quantum entanglement between the two subsystems, can be measured by a 
well-known quantity: the quantum linear entropy |l(1 | 

S q = l-Tn05 2 ), (2) 

where Tr, denotes a trace over the i-th subsystem, and p = Tr 2 (\ip(t))(ip(t)\) is the reduced density operator for the 
first subsystem. An increase in S q suggests an increase of 1/(1 — S q ), which gives the number of orthogonal quantum 
states that are incoherently populated if the second subsystem is ignored. Below we choose q,p as the "bath" variables 
and P, Q as the system variables. 

Let p c (Q, P, q,P, t) denote the phase space distribution function evolved classically, and pw(Q, P, <7,P, t) denote the 
quantum (Wigner) phase space distribution function. Their time evolution equations are given by 

7^ = {H,Pch (3) 

= {H,p w } M , (4) 

where {•} denotes classical Poisson bracket and {-}m denotes quantum Moyal bracket We define classical and 
quantum reduced distribution functions as 

p c (Q,P,t) = J p c (Q,P,q,p,i) dq dp, (5) 

pw(Q,P,t) = J pw(Q,P,q,P,t) dq dp. (6) 

Since 

S q (t) = 1 - 2nh J p 2 w (Q, P, t) dQ dP, (7) 

where h is the effective Planck constant, we can define a classical analog [denoted S c (t)] to S q (t) by replacing pw 
with p c . That is, 

S c (t) = l-2irh J pl(Q, P, t) dQ dP. (8) 

The main focus here is to compare S c (t) with S q (t), i.e., the classical vs. quantum evolution of the intrinsic decoherence 
dynamics, as measured by the classical vs. quantum entropy. 

Perturbative treatments have proved to very useful in understanding decoherence dynamics |sl IT2I fill. Il4| . Here, to 
analytically examine classical vs. quantum intrinsic decoherence dynamics at early times, we apply the perturbative 
approach developed in our previous work 8] to the case of intrinsic decoherence dynamics. Specifically consider a 
second-order perturbative expansion with respect to the time variable t for both S q and S Cl i.e., 

t t 2 

S c {t) = S c (0) + — + — + ... , 

T c,l T c2 

S q (t) = S q (0) + J- + -*- + •■• . (9) 

T q.! T q,2 
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Then, from the classical and quantum dynamics of the entire system one obtains 

-L = -47Tft J p c (Q,P,0) J {H,p c (Q,P,q,p,0)} dqdpdQ dP, 



(10) 



-L = -ink J p w (Q,P,0) J {H,p w (Q,P,q,p,0)} M dq dp dQ dP, 



(11) 



and 



1 

5 



1 

r q,2 



-2nh J p c (Q,P,0) J {H,{H,p c (Q,P,q,p,0)}} dqdpdQ dP 
-2nh J J {H,p c (Q,P,q,p,0)} dq dp 



dQ dP, 



-2nh J p w (Q,P,0) J {H,{H,p w (Q,P,q,p,0)} M } M dq dp dQ dP 
2lTh I I { H,Pw (® ,P,q,p,() ^ M dqdp 



dQ dP. 



(12) 



(13) 



Further, using the definitions of the classical Poisson and quantum Moyal brackets and assuming that initial classical 
and quantum distribution functions are identical and separable, i.e., 



we have 



p c (Q,P,q,p,0) = p w (Q,P,q,p,0) = ft(Q,P)f%(q,p), 



1 1 

T C ,1 7"q,l 



0. 



(14) 



(15) 



and 



= 2irh 



1 



'c,2 



1 

<2 



dp\{Q,P) 



-i 2 



dP 



C(0, 0) dQ dP, 



2nh 



h^h>0 



[n/(2i)} 



(2h+2h) 



(2h + l)!(2/ 2 + 1)! 



d( 2h+r >fi(Q,P) d^+^p\{Q,P) 



Qp(2h + l) 

where C(h, I2) is a correlation function given by 



Qp(2l 2 +1) 



C(h,l 2 ) dQ dP, 



C(h,l 2 ) 



/a (j"i+i)y(Q ig ) d^ h +^V(Q,g) 
gQ(2i!+i) dQ(2*2+i) 



d(^ + 1 )V{Q,q)\ /d^+^V(Q,q) 



9Q(2h+i) /.„ \ 9Q(2/ 2 +i) , ;/ , 



(16) 



(17) 



(18) 



Here (-)^o denotes the ensemble average over the zero-time "bath distribution function" p^q^p). It is worth empha- 
sizing that in our derivations we have used the same initial state for the classical and quantum dynamics. 

Equation (|15fl shows that the zero first-order linear entropy increase rate, i.e., l/r q i=0, has a strict classical 
analog. Further, Eq. i|16fl indicates that classical Liouville dynamics also predicts a second-order entropy production 
rate 1 /t% 2 that is the analog of the second-order quantum decoherence rate 1/T q 2 . Thus, we can identify two categories 
of early-time intrinsic decoherence dynamics: classical if r Ci 2 ~ T q .2, and nonclassical if r q .2 appreciably differs from 

T c ,2- 
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To simplify Eqs. (|16(l and l|17|) we introduce the Fourier transform [denoted F(Qi,Q 2 )} of f>1(Q,P), i.e., 

~iAQP~ 



F(Q U Q 2 ) 



ft(Q,P)exp 



dP, 



where Q = (Qi + Q 2 )/2 and AQ = Qi — Q 2 - We then obtain 

-5- = i / I^(Qi,Q 2 )| 2 aq 2 c(o,o) dQi dg 2 

T c : 2 ft J 

and 

^g(2/!+2; 2 +2) J 



Ill 

^1 ~ <^ + ft 5 



/ ito,<mi' £ <^ 



(2Zi + 1)!(2Z 2 + 1)! 2( 2 ^+ 2 '2> 



C{h,l 2 ) dQi dQ 2 . 



(19) 



(20) 



(21) 



Equations (|20|l and (|21|) are general results. For the special case of Vi2(Q, q) = f(Q)g(q), Eqs. I|2U|I and (|21|l can be 
rewritten in a simple and more enlightening form: 



1 (g 2 (d))p° (9(d))l« 



'c,2 



|f(Qi,q 2 )| 2 aq 2 



d/(Q) 



and 



1 (5 2 W)pO - (g(q))lo 



'q,2 



|F(o l7 g2)| 2 AQ 2 



A/(Q) 
AQ 



(22) 



(23) 



where Af(Q)/AQ is the finite-difference function 



A/(Q) _ f(Q + AQ/2)-f(Q-AQ/2) 



AQ 



AQ 



/(Qi)-/(Q 2 ) 
Q1-Q2 



(24) 



As a result: (1) If f(Q) depends only linearly or quadratically upon the coupling coordinate Q, a common approxi- 
mation, then (l/r 2 2 — 1/t 2 2 ) = for any initial state. That is, in this case there exists perfect QCC in early-time 
dynamics of intrinsic decoherence, regardless of h, and irrespective of the potentials Vi(Q) and V^q). (2) Even in the 
case of highly nonlinear /(Q), as long as F(Q 1: Q 2 ) decays fast enough with \Qi — Q 2 \ such that A/ / AQ « df/dQ, 
QCC would still be excellent. The smaller the H, the more rigorous is this requirement. (3) If Af/AQ differs sig- 
nificantly from df/dQ over the Q-coordinate scale of the initial state, quantum entropy production can be totally 
unrelated to classical entropy production. Such cases of poor QCC are of fundamental interest, but are not the focus 
of this paper. 

The second-order perturbative treatment is most reliable for early-time dynamics and for relatively weak decoher- 
ence. The above results are particularly significant for studies on the control of intrinsic decoherence, where early-time 
dynamics of weak decoherence is important. In these circumstances it is useful to understand the extent to which 
(quantum) intrinsic decoherence is equivalent to classical entropy production, i.e. to increasing S c (t). In particular, 
if there exists good correspondence between classical and quantum decoherence dynamics, then the essence of deco- 
herence control is equivalent to the suppression of classical entropy production, and various classical tools may be 
considered to achieve decoherence control. If not, then fully quantum tools are required. 

The above perturbation results clearly demonstrate that quantum dynamics of intrinsic decoherence has a direct 
analog in classical Liouville dynamics. This rather intriguing result motivates us to computationally examine QCC 
in the dynamics of intrinsic decoherence over all time scales. 



III. COMPUTATIONAL RESULTS: TWO SAMPLE CASES 



To computationally examine QCC in the dynamics of intrinsic decoherence, we consider coupled-oscillator model 
systems with smooth Hamiltonians. In all the model systems studied below, we choose 



V 1 (Q) + V 2 ( q )^^-(Q 4 + q 4 ) 



(25) 
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where j3 = 0.01. Since Vi(Q) and V^q) have no simple harmonic terms, any observed agreement between classical and 
quantum behavior cannot be attributed to the similarity between classical and quantum harmonic oscillator dynamics. 
If the coupling potential Vi2(Q, q) is quadratic in both Q and q, i.e., VWQ, q) = aQ 2 q 2 /2, then the resultant coupled- 
oscillator system is the well-known quartic oscillator model [lj, iWL i1T\ . Because this model is well-studied and can 
display strongly chaotic (e.g., a — 1.0, (5 — 0.01) or integrable (e.g., a — 0.03, (3 — 0.01) dynamics, it is used in Sec. 
IIVI as an ideal model to study QCC in intrinsic decoherence dynamics for both integrable and chaotic cases. 

Our perturbation theory approach predicts good classical-quantum agreement at short times for some potentials 
and initial conditions and poor agreement for others. We examine both these cases computationally. 

It suffices to consider one case of poor agreement, since poor QCC at early times invariably translates to similar 
behavior at later times. Consider then Vi2(Q, q) to be some highly nonlinear potential. Computations of the quantum 
dynamics and thus the time dependence of S q (t) are straightforward ^(|. S c (t) is computed directly using Monte- 
Carlo simulations with an importance sampling technique (where the Monte-Carlo simulations are based upon Eq. 
(|27H below). From the analytical results above we see that Vi2(Q, q) and the scale of the initial state play decisive roles 
in QCC in early-time intrinsic decoherence dynamics. In particular, we expect poor QCC if Vi2(Q,q) — f(Q)g(q) 
differs significantly from a linear or quadratic function of Q such that A//AQ differs significantly from (df/dQ) 
over the Q-coordinate scale (i.e., the support) of the initial state. To confirm this computationally we consider 
f(Q)g(q) = sin 2 (10(2)q 2 , with the initial distribution functions of the two subsystems given by 



P°(Q,P) 



— exp 



(Q-Q0) 2 (P-Pof 



7i h 



exp 



(q - qo) 2 (p - po) 2 



2a* 



2al 



(26) 



Here the dimensionless effe ctive Planck const ant i s chosen to be % 
and ctq/25 = 25 • <j p = a q = a v = y/h/2, Q = 0.5, P a = 



= 0.005 throughout, except for one case in Sec. 
0.5, q = 0, with H(Q ,P ,q ,po) = 0.24. Note 



that pi(Q,P) is strongly squeezed in P and that this initial distribution function is considerably delocalized in Q. 
Further, since |/(Q)| = | sin 2 (10Q)| < 1.0, |AQ| • \df/dQ\ can be much larger than |A/(Q)|. Thus, for this case the 
perturbation result predicts that at early times there can be substantial classical entropy production with insignificant 
quantum decoherence. As shown in Fig. ^ this is nicely confirmed by the numerical results of S q (t) and S c (t). In 
particular, Fig. shows that, at t = 1.0 S c {t) (discrete points) is ~ 0.9 while S q (t) (solid line) is still less than 0.2. 
Evidently, QCC in this case is indeed very poor from the very beginning. 




FIG. 1: A comparison between Sq(t) (dashed line) and S c (t) (discrete circular points) in the first sample case. The coupling 
potential is highly nonlinear such that at early times classical entropy production is much faster than quantum entropy 
production. See the text for details. All variables are in dimensionless units. 



There remains then the important question of the quantitative degree of QCC in circumstances where our per- 
turbative analysis predicts good short-time QCC. In particular, it is important to investigate whether or not good 
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QCC predicted perturbatively remains for a considerable amount of time. If so, then the perturbative treatment 
provides a useful guide to our understanding of QCC in intrinsic decoherence dynamics. If not, then our perturbative 
results make sense only for extremely weak decoherence. Dramatically, our computational studies strongly support 
our analytical perturbation results, even in the presence of significant decoherence. For example, consider the case, 
where the parameters for the initial state are the same as in the previous case (therefore the initial state is also much 
delocalized) , but the coupling potential is given by Vi2(Q, q) = Q 2 sin 2 (g). This coupling potential is highly nonlinear 
in q but still quadratic in Q. In accord with the second-order perturbation results, such a coupling potential should 
still give rise to good early-time QCC in the intrinsic decoherence dynamics of the first subsystem. This is confirmed 
by the quantitative comparison between S q (t) and S c (t) shown in Fig. More importantly, Fig. |2] shows that 
outstanding QCC remains even when both S q (t) and S c (t) have increased to close to their saturation value of unity. 
Numerous other computational results (not shown) are consistent with the two cases shown here. 




FIG. 2: A comparison between <Sq(i) (dashed line) and S c (t) (discrete circular points) in the second sample case. The coupling 
potential is highly nonlinear in terms of the position of the second subsystem, but is quadratic in terms of the position of 
the first subsystem, resulting in excellent quantum-classical correspondence in intrinsic decoherence dynamics even though the 
initial distribution function of the first subsystem is considerably delocalized. See the text for details. All variables are in 
dimensionless units. 

These results show the usefulness of the second-order perturbation theory in understanding QCC in intrinsic deco- 
herence dynamics emanating from squeezed initial states. Also of interest is intrinsic decoherence dynamics associated 
with sufficiently localized initial states, which, in accord with the previous perturbation results, should display excel- 
lent early-time QCC for any coupling potential V12 (Q,q). We now computationally examine QCC at later times for 
localized states as initial conditions and explain the results in terms of a simple classical theory of intrinsic decoherence 
dynamics. 



IV. LOCALIZED INITIAL STATES 



A. Simple Classical Approach 

Below we show that S q (t) and S c (t) are often in excellent agreement, over large time scales, for initially localized 
states, significantly extending the perturbation theory result. In doing so we compare the quantum S q (t) with full 
classical mechanics as well as with a simple classical theory derived in this section. The latter provides further insight 
into the origins of increasing S c (t). 
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To derive the simplified classical result we first use Liouville's theorem to reexpress S c (t) as 

S c (t) = l-2izhj J Pc(Q(t),P(t),q(t),p(t),t) 

Xp c (Q(t),P(t),q',p',t) dq' dp' dQ(t) dP{t) dq{t) dp{t) 
= 1-2%% J J p c (Q,P,q,p,0)p c (Q(t),P(t),q',p',t) dq' dp' dQ dP dq dp 

I J p c (Q,P,q,P,0)Pc(Q",P",q",p",0) dq' dp' dQ dP dq dp, 



= 1-2tt7L 



(27) 



where (Q(t), P(t),q(t),p(t)) is the phase space location of the trajectory emanating from (Q,P,q,p) at t = 0, and 
(Q" , P" , q" ,p") is the phase space location of the trajectory at time zero if the classical trajectory is propagated 
backwards from (Q(t), P(t),q' ,p'). Because the initial state p c (Q,P,q,p,0) is assumed highly localized in phase 
space, (Q" , P" ,q" ,p") must be very close to (Q,P,q,p) in order for the term p c (Q, P,q,p,0)p c (Q" , P" , q" ,p" ,0) in 
Eq. (|27[) to be appreciable and thus to contribute to S c (t). Hence a convenient approximation can be made: we 
assume that, at time t, only those backward trajectories near (Q(t), P(i), q(t),p(t)) need be taken into account. This 
means that we treat Q" — Q, P" — P, q" — q, p" — p, Sq' = \q' — q(t)], and 5p' = [p' — p(t)] as sufficiently small such 
that 



P" 



Q 
p 

q- 
p- 



-M 13 (t)6q' 
-M 23 (t)Sq' 
M 33 (t)Sq' 4 
M 43 (t)VH 



- M u (t)5p\ 
-M 2i (t)S P ', 
M M (t)5p', 
M M (t)Sp', 



(28) 



where My = 1,2,3,4) is the stability matrix associated with the backward trajectories emanating from 

(Q(t),P(t),q(t),p(t)): 

d(Q,P,q, P ) 



y d(Q(t),P(t),q(t),p(t))- [ ' 

Although this approximation should be less reliable for chaotic systems, we demonstrate below that, it is, nonetheless, 
computationally useful in both integrable and chaotic cases. 

For simplicity we consider below a specific example in which p\(Q,P) and p^iq^p) are symmetric Gaussian states 
given by Eq. lf^|) with aq = op =a q = a p = a = \Jfij2. Other distributions can be considered in an analogous 
fashion. Substituting Eqs. I|26() and (|28|l into Eq. I|27|) and evaluating the integrals, we obtain 



1 / exp 



t<-X 2 +V 1 Y J -2UVZ 



2a 2 (X 2 Y 2 -Z 2 ) 



VX 2 Y 2 - Z 2 



(30) 



where 



X = 


M* 3 + Ml 3 + Ml 3 + Mi 3 , 






Y = 


+ M| 4 + Mf 4 + M| 4 , 






Z = 


M 13 M 14 + M23M24 + M33M34 + M43M44, 






u = 


(Q - Q Q )M 14 + (P - P )M 24 + (q- q )M 34 + 


(P 


- pa)M 44: 


V = 


(Q - Q )Af 13 + (P - P )M 23 + {q - q )M 33 + 


(P 


- p a )M i3 




p' c (Q, P, q, p) = 4n 2 h 2 [p?(Q, P)] 2 [p° (q, p)} 2 







(31) 

and 

(32) 

Equations (|30|l and l|31|l indicate that the classical dynamics of intrinsic decoherence is closely related to the classical 
stability matrix elements averaged over a rescaled initial distribution function. This interestin g co nnection provides 
insight into a variety of interesting aspects of quantum intrinsic decoherence dynamics 0, Il9l l2fl l2ll 122^ . For 
example, for chaotic dynamics in which classical trajectories are highly unstable and therefore in which |My | increases 
rapidly, S c (t) should increase much faster than for the case of integrable dynamics. This observation can, with the 
assumption that there is fairly good QCC in intrinsic decoherence dynamics, directly explain previous results on 
quantum signatures of classical chaos in the dynamics of quantum entanglement |18|. Further, because Eqs. I|30|l and 
(|31|l are expressed in terms of classical stability matrices, characteristics of the time dependence of S c (t) and therefore 
of Sq(t) can be easily related to the time and space fluctuations in the instability of classical trajectories. 
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B. Computational Results: Localized Initial States 




time 



FIG. 3: A comparison between Sq(t) (dashed line) and the approximate S c (t) (solid line) calculated from Eq. I3UI for the 
quartic oscillator model in the case of integrable dynamics (a = 0.03, f3 — 0.01). The initial state is given by Eq. 12611 . with 
<yp = o-q = <7p = a q — h = 0.005, Qo = 0.4, Po = 0.5, qo = 0.6, and H(Qo, Po, qo, Po) = 0.24. Full classical results based 

upon Eq. (1271 are represented by discrete circular points. All variables are in dimensionless units. 




time 



FIG. 4: Same as in Fig. ^except for strongly chaotic dynamics (a = 1.0, j3 = 0.01). 

Consider then QCC over large time scales for localized initial states for both integrable and chaotic cases. To 
do so we examine the quartic oscillator model as well as results where the coupling potential is replaced by the 
nonlinear potential V^iQ, q) = Q.5Q 2 q 2 + Q A q 2 ■ The initial states are chosen to be localized initial states, and both 
the full classical dynamics and the approximate time dependence of S c (t) in Eq. I|3U|) are compared to the quantum 
result. Specifically we realize the ensemble average in Eq. I|30(l by Monte-Carlo simulations, using only 2 x 10 4 
sampling classical trajectories from which the stability matrix elements Mij are evaluated. The initial Gaussian 
states are chosen symmetric with a — y/%/2 — 0.05 and are sufficiently localized so that Eq. l(2*Hjl should be a valid 
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approximation. Full classical results for S c (t) (represented again by discrete circular points), which are much more 
demanding computationally, are also provided below. 

Figures Eland 01 compare results for S c (t) obtained from Eq. H30|l and from exact classical results with S q (t) (dashed 
line) for integrable and chaotic dynamics in the quartic oscillator model, respectively. A number of observations are 
in order. First, it is clear that in both cases the approximate S c (t) are in excellent agreement with the full classical 
results, confirming the utility of the simple model [Eq. Second, both S c (t) and S q (t) are seen, in the chaotic 

case, to relax faster towards 1.0 than they do in the integrable case. Third, the oscillation amplitudes of S q (t) in 
the chaotic case are much smaller than that in the integrable case. Hence, the fast relaxation and small-amplitude 
oscillations of S q (t) shown in Fig. 01 may be regarded as fingerprints of the underlying classical chaos. Finally, and 
most importantly, the entire time dependence, including oscillations in Fig. |3| and Fig. 01 of S q (t) are beautifully 
captured by both the exact and the approximate S c (t). It should also be noted that the QCC time scale shown in Fig. 
0]is appreciably longer than is the QCC break time tb ~ 5.0 for the same Ti, obtained by quantitatively comparing the 
structure of the classical and quantum distribution functions |2^. This can be understood by the fact that S c (t) [or 
S q (t)] describes the reduced distribution functions p c (Q, P, t) [or pw(Q, P, t)} which is insensitive to the fine structure 
of p c (Q,P,q,p,t) [or p w (Q,P,q,P,t)}. 

Calculations for many other initial states confirm that the QCC results shown in Figs. |21and0]are typical, indicating 
that (a) QCC is essentially exact over large time scales and (b) the simple classical theory of intrinsic decoherence 
dynamics introduced above provides a useful approximation to the exact results. 
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FIG. 5: Same as in Fig. except for strongly chaotic dynamics (a 
with Q = P = 0, g = 0.6 and H{Q , P , q ,Po) = 0.24. 



1.0,/? = 0.01), and for a special initial Gaussian state, 



Figure shows one case, however, where the approximate S c (t) and exact classical or quantum results differ 
quantitatively. Here the system is still the quartic oscillator model with a = 1.0, (3 = 0.01, but with an initial state 
of special type. In particular, both the initial average position Qo and the initial average momentum Pq are set to 
zero. Initial states of this type are called channel states 01 > an d effectively give rise to very weak coupling between 
the two subsystems over a considerably large time scale. Indeed, since at short times df/dQ » Af(Q)/AQ as 
for Q ps Qo = 0, one obtains from Eqs. (|22|l and (|23l) that 1/t^ 2 = 1/ T lq ~ 0- Therefore the early-time intrinsic 
decoherence rate should be small, as seen in Fig. 03 although the underlying classical dynamics is strongly chaotic. 
For this reason one expects that the dynamical behavior of S c (t) and S q (t) should differ from previous cases. As 
shown in Fig. in this case both S q (t) and S c (t) increase in a step-wise fashion, distinctly different from that in Figs. 
and 01 Agreement between them remains excellent. However, the approximate S c (t) misses some of the important 
structure. 

In Fig. EJwe show the QCC result for a simple variant of the quartic oscillator model, i.e., Vi2(Q, q) — 0.5Q 2 q 2 +Q 4 q 2 , 
a coupling potential that is neither linear nor quadratic. As seen in Fig. El even with such nonlinear coupling S c (t) and 
S q (t) are in excellent agreement over large time scales. This emphasizes the fact that the good QCC results observed 
in the quartic oscillator model are not due to the fact that the coupling potential therein is quadratic. Hence, we 
conclude that for initially localized states, our simple classical theory of intrinsic decoherence dynamics [see Eq. (|30l) ] 
is generally useful in describing intrinsic decoherence dynamics in smooth Hamiltonian systems. 
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FIG. 6: Same as in Fig. EJexcept for a modified quartic oscillator model in which V12 (Q,q) = 0.5Q 2 q 2 + Q 4 q 2 



V. DISCUSSION AND SUMMARY 



Quantum entanglement between individual subsystems has no classical analog. Nevertheless, as shown in this work, 
the quantum dynamics of quantum entanglement, as manifest in the quantum dynamics of intrinsic dccohcrence, does 
have a classical analog in classical Liouville dynamics describing classical correlations between classical subensem- 
bles. Hence it is useful to isolate the conditions under which there is good QCC in intrinsic decoherence dynamics. 
This is done analytically, for early-time dynamics and for weak decoherence, by a second-order perturbative theory. 
Interestingly, as demonstrated by our computational studies in Sec. IIIII and Sec. IIVI the physical picture of QCC 
afforded by the perturbative treatment can be still very useful even when the time scale under investigation is rela- 
tively long and the degree of intrinsic decoherence is significant. In particular, under the circumstances where there is 
good early-time QCC, classical Liouville dynamics can provide a simple means of understanding different aspects of 
intrinsic decoherence dynamics, for relatively large time scales and for both integrable and chaotic dynamics. Further, 
we have derived an approximate but very simple classical theory of linear entropy production of intrinsic decoherence 
dynamics associated with localized initial states, and shown that the rate of entropy production is closely related to 
the stability properties of classical trajectories. 

Clearly, the linear entropy is just one of many possible representation-independent measures of intrinsic decoherence, 
and S q (t) w S c (t) does not mean that the quantum dynamics is equivalent to the corresponding classical Liouville 
dynamics. For example, if the saturation value of the linear entropy in the long time limit is of particular interest, 
then the measures 1/[1 — S c (t)] and 1/[1 — 5 q (£)] (which gives the number of orthogonal states that are incoherently 
populated) should be more useful in describing QCC. Indeed, our results in Figs. 2]and[|)]suggest that as time increases 
one has 1/[1 — S c (t)] >> 1/[1 — S'q(i)]. This is consistent with our previous observation [8| that, decoherence can 
dramatically improve QCC, but even strong decoherence does not necessarily suffice to ensure that quantum entropy 
production is the same as classical entropy production. 

It should also be pointed out that the model quantum systems studied in this paper are still far from the semiclassical 
regime. This is indicated, in the chaotic case of the quartic oscillator model for example, by the fact that the QCC 
break time is relatively short compared to the time scale that we examined. Correspondence will worsen quantitatively 
with increasing h, although, as discussed above, h is far from the only factor influencing the quality of the QCC. 
However, we note that for localized initial states the qualitative features of the time dependence of classical and 
quantum linear entropy may remain similar to one another with much larger effective Planck constants. For example, 
Fig. [7| displays fairly good QCC between S q (t) and S c (t), in the chaotic case of the quartic oscillator model, with 
h = 0.05 and with an initial symmetric Gaussian state. 

A number of interesting extensions of this work are under consideration. First, it seems straightforward but 
necessary to consider cases in which the coupling potential depends upon both position and momentum. Second, 
we propose to further investigate the role of the dynamics of the subsystems in addition to that of the coupling 
potential (e.g., the dynamics in coupled Morse oscillator systems). Third, it is interesting to study QCC in intrinsic 
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FIG. 7: A comparison between S q (t) (dashed line) and S c (t) (discrete circular points) for strongly chaotic dynamics of the quartic 
oscillaor model (a = 1.0, (3 — 0.01) and for h = 0.05. The initial state is given by Eq. 1261 . with op — oq = a v — a q — 
Qo = 0.4, Po = 0.5, qo = 0.6, and H(Qq, Po, qo,Po) = 0.24. All variables are in dimensionless units. 

decoherence dynamics in terms of the decay of off-diagonal density matrix elements. Such studies are ongoing, with 
preliminary studies [^[indicating that comparing the time dependence of off-diagonal density matrix elements to its 
direct classical analog 8] will provide deeper insights into QCC in intrinsic decoherence dynamics. 

To summarize, we have shown that classical dynamics can be very useful in describing intrinsic decoherence dynamics 
in smooth Hamiltonian systems. In particular, we have identified conditions under which excellent quantum-classical 
correspondence in the early-time dynamics of intrinsic decoherence is possible via a second-order perturbative treat- 
ment, have presented a simple classical theory of intrinsic decoherence dynamics emanating from localized initial states, 
and have provided supporting computational results. The hope is that by extending this study to high-dimensional 
Hamiltonian systems, we may use purely classical approaches to describe (at least qualitatively) the dynamics of 
quantum entanglement or intrinsic decoherence in polyatomic molecular systems. 
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